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Abstract 

The instability of nuclear matter due to particle-hole excitation modes has been studied in the 
frame-work of several relativistic mean field (RMF) models. It is found that both the longitudinal 
and the transversal modes depend sensitively on the parameter sets used. The important impact 
of the vector and vector-scalar nonlinear terms on the stability of both modes is demonstrated. 
Our finding corroborates the result of previous studies, namely that certain RMF models cannot 
be used in high density applications. However, we show that for certain parameter sets of RMF 
models this shortcoming can be alleviated by adding these nonlinear terms. 

PACS numbers: 21.30.Fe, 21.65. +f, 21.60.-n 
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Remarkable progress has been made in constraining the nuclear equation of state (EOS) 
from astrophysics and heavy ion reactions, i.e., we know at a confident level that the EOS 
should be soft at moderate densities but relatively stiff at high densities (for recent paper 
see, e.g., Ref {l']). However, the fast development in the unstable nuclear beam facility ^ 
will also reveal a lot of unexpected phenomena in unstable nuclei far from the stability line 
region in the near future. Therefore, a unified model describing simultaneously finite nuclei 
and matter properties at high densities with high degree of accuracy is not only challenging 
but also mandatory. 

Relativistic mean field (RMF) models have been quite successful in providing a micro- 
scopic description of many ground states properties ranging from medium, heavy, up to 
super- heavy nuclei (for a review see, e.g., Ref {s]). The standard ansatz (S-RMF) uses a, 
uj and p mesons as degrees of freedom with additional cubic and quartic nonlinearities of 
a meson to describe the interaction. The corresponding parameter sets of this model are 
known as, e.g., NL-Z {4] and NL3 5]. The simplest extension of the S-RMF is achieved 
by introducing a quartic nonlinearity of the u meson in the Lagrangian (V-RMF). The pa- 
rameter sets TMl and PKl belong to this model version. Other extensions (E-RMF) 
of the S-RMF parameterization are known as Gl and G2 [8]. The corresponding model is 
derived from the effective field theory which allows for possible scalar- vector coupling terms 
up to fourth order to be present in the Lagrangian. Other interesting properties of the RMF 
models come from the fact that the relativistic nature of the models and the properties 
of nuclear matter at saturation are fulfilled while in the extrapolation to higher densities 
the appearance of acausal behavior (the speed of sound exceeds the speed of light) can be 
alleviated [Oj. 

Surprisingly, though, little attention has been given so far to check the matter instability 
by means of particle-hole excitations with frequency ^0=0 in the transversal and longitudinal 
modes, based on these models at high densities. It is understood that this analysis is one of 
the tools to check the reliability of the models at high densities. Efforts in this direction have 
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been devoted some times ago for the linear model case pXJ, , whereas the latest progress 



can be found in Ref. [12] , in which the S-RMF model was used to investigate the instability 
of symmetric nuclear matter (SNM). The authors of Ref. [l2| found that the onset of the 
instability depends sensitively on the parameterization of the model. Parameterizations 
with a low nucleon effective mass {m*/m ~ 0.5-0.6), which can accurately predict nuclear 
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ground state properties, are critical in the longitudinal mode. The S-RMF model produces 
also a wide instability regime in the transversal mode, even for the parameterization with 
high nucleon effective mass or by using a more general functional form of nonlinear a- 
meson self-interactions. In addition, it is known that at the critical density, which is larger 
than saturation density, the a meson self-interaction in the S-RMF model becomes unstable 
because the square of the effective o meson mass becomes negative and, as a consequence, 
an additional instability regime in longitudinal mode appears (see region II in the lower-left 
panel of Fig. [1]). This instability is not present in the linear RMF model. To overcome the 
mentioned shortcomings, the authors of Ref. 12] suggested to add nonlinear terms which 
contain not only functions of the a meson field, but also of the uj meson field. 

The analysis of the instability in symmetric nuclear matter (SNM) is a first step towards a 
more comprehensive analysis taking into account matter with multi-component constituents, 
where the latter plays a key role in our understanding of some stellar matter problems. 
Due to the simplicity of SNM, many essential physical aspects can be understood lucidly. 
Therefore, we will revisit the instability problem of the RMF models at high-density SNM 
but now the role of vector and vector-scalar coupling nonlinearities are taken into account 
and the effects in the unstable regimes with respect to longitudinal and transversal modes 
are investigated by using some selected parameterizations of S-RMF, V-RMF and E-RMF 
models, where the S-RMF is used as the benchmark. 

To determine the instabilities of the RMF models, we start from the energy density e in 
SNM which takes the following form, 

- \d,V^ - d^aV^ - \d,a^Vl (1) 

where for the NL-Z and NL3 parameter sets (S-RMF) the last three terms vanish, and for 
the TMl and PKl parameter sets (V-RMF) only the terms proportional to d2 and d^ vanish, 
while for Gl and G2 parameter sets (E-RMF) all parameters are utilized, eunear is a function 
of the kinetic terms of the nucleons, a and uo masses as well as interactions terms of the 



nucleons. Using a similar procedure to the one used in Refs. 10|, ll3| , we can calculate the 



transversal (er) and longitudinal (e^) dielectric functions of the RMF models. For SNM, 
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they have simple forms, i.e., 



6i = 1 + 2dsUs - 2d^Uv + M^sv'^sv, (2) 



with 



Uv = Hoo - 2dsnl, + 2dsUsUoo 
Usv = Hm + 2d^sv^l + 2dsv'ns'noo, (3) 

where 11^, Us, ^oo, and 11^/ are the transversal, scalar, longitudinal and scalar-vector cou- 
pling polarizations of proton or neutron with go = 0. In SNM, the values of each polarization 
for proton and neutron are equal. The explicit forms of these polarizations are given in 
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Refs. [10|, 113| . The longitudinal scalar meson propagator is given by 

ds = (4) 

while the vector meson longitudinal and transversal propagators are 



dy 



q2 _|_ ^* 



- o , 1,, 2' (^) 



and the scalar-vector coupling propagator takes the form 



where m*J= d'^e/da'^, ml^= -d^e/dVo^ and ^^^=-d'^e/ dadVQ. It is clear that the NL3, 
NL-Z, TMl and PKl parameterizations have vanishing A^tj and only S-RMF parameteri- 
zations have a constant effective omega-meson mass, which is equal to its bare mass. The 
unstable regimes are determined from e/,, < 0. The corresponding results for some se- 
lected parameter sets are shown in Fig. [H where the left panels exhibit the longitudinal 
modes and the right panels display the transversal ones. 

The NL3 and NL-Z models have two instability regimes (I and II). In regime II, NL3 has 
a deeper valley compared to NL-Z and for each parameter set, the valley in regime II is 
deeper than regime I. Regime II of the NL-Z parameter set starts to appear at pc ~ 2.5 
while for NL3 it starts at pc ~ 6.5 po) both with qc < 2.5 fm. If we use an artificial parameter 




FIG. 1: (Color online) Longitudinal (left panels) and transversal (right panels) modes of the 
instability due to the particle-hole excitation with qq = for the SNM and some selected RMF 
parameterizations. Different depths of the dielectric functions are shown by different colors (grey 
scales) . 

set with a large nucleoli effective mass {m*/m > 0.7 at saturation), the unstable regime II 
disappears. On the other hand, for the V-RMF and E-RMF models with m* /m ~ 0.5-0.6 
(the effective nucleon mass range for which the correct spin-orbit splittings is reproduced) 
at saturation, the regime II does not also exist. For NL-Z, regime I appears quite early, i.e. 
at pc ~ 2.5 po with a relatively small qc (~ 2 fm), while for NL3 it begins at pc ~ 5 po with 
gc ~ 3 fm. 

The vector and vector-scalar coupling in the nonlinear terms of the V-RMF and E-RMF 
models lead to somewhat narrower unstable regimes compared to the standard one. For 
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FIG. 2: (Color online) Longitudinal scalar (ds), vector (dy) and vector-scalar coupling propagators 
(dy^)(left panels) and the corresponding contributions to the longitudinal dielectric function (right 
panels). Here g = 3 fm is used. 

TMl, it appears at pc ~ 10 po with gc ~ 5 fm, for PKl it starts at pc ~ 5 po with qc ~ 3.5 
fm, while for G2 it starts at ~ 22 po with ~ 6.5 fm and for Gl it appears at p^ ~ 3 
Po and gc ~ 3 fm. This means that in the V-RMF and E-RMF models there are parameter 
sets, e.g. TMl and G2, for which their instabilities in the longitudinal mode can be shifted 
into the regimes which are physically not too important (large qc and Pc) without losing 
their accurate predictions in ground-state properties of finite nuclei. 

Contributions of the longitudinal-scalar {ds), vector (dy) and vector-scalar {dyg) prop- 
agators in the longitudinal dielectric function of Eq. ([2]) with g = 3 fm (below q^ of 
V-RMF and E-RMF models) for all parameter sets used can be seen in Fig. [2l The den- 
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FIG. 3: (Color online) Same as Fig. [2] but with q = 0.5 fm. Different with other parameter sets 
which have small regular contributions, contribution of the NL3 and NL-Z parameter sets at certain 
density scalar propagator exhibit discontinuities (lower left) which leads to the appearance of the 
large positive dylly (middle right) and negative dg^s (lower right) above this density. 

sity dependence of the uj meson propagator in the V-RMF or E-RMF models produces a 
negative vector contribution {dyllv) for all densities. In the case of the Gl parameter set, 
the vector-scalar coupling {dgylisv) correction enhances the stability, but this contribution 
is smaller than Gl albeit with different behavior for G2. Overall, these yield a sufficient 
suppression to the negative contribution from the scalar one [dsTls) at low density regime. 
As a consequence, a positive is produced. In contrary, for the S-RMF model, a constant 
UJ meson propagator allows relatively large positive vector and negative scalar contributions. 
This opens the possibility that becomes negative at low densities. 

The same contributions but with q = 0.5 fm can be seen in Fig. [31 The square effective a 
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FIG. 4: (Color online) Transversal propagators with g = 5 fm (lower-left panel) and q = 0.5 fm 
(upper- left panel) followed by their transversal dielectric functions (right panels). 

meson mass of the NL-Z and NL3 parameter sets becomes negative after reaching a certain 
critical density (pc) due to particular parameter values of their scalar nonlinear terms, which 
leads to the appearance of a discontinuity in a meson propagator at that point. For NL3, the 
effect is more dramatic but the regime is shifted to a higher critical density compared to the 
NL-Z parameter set. Furthermore, large positive vector and negative scalar contributions 
for p > Pc appear in this model. Therefore, in this density range, the longitudinal mode 
becomes unstable, even for very small momentum response {q ^ 0). This is the reason that 
for both parameter sets the unstable regime II exists. 

The S-RMF model has a broad instability for this mode. It does not depend on the pa- 
rameterization used, except for low momentum response {q < 1 fm). The interesting finding 
here is that for TMl and Gl parameter sets the unstable regimes absolutely disappear. The 
particular parameter values of the vector and vector-scalar nonlinear terms are responsible 
for this case. While for the PKl and G2 parameter sets relatively narrow unstable regime 
are produced, the regimes are also shifted to a relatively high density and large momentum 
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FIG. 5: (Color online) Unstable longitudinal regimes I for asymmetric nuclear matter using NL3 
parameter set with proton fraction equals 0.4 in the left panel and 0.22 in the right panel. 

( Pc ~ 5 po with gc ~ 3 fm for PKl and Pc ~ 8 po with ~ 4 fm for G2). This behavior 
never occurs for the S-RMF model parameterizations. The corresponding results for some 
selected parameter sets are shown in the right panels of Fig. [1] The role of the density- 
dependent effective-transversal omega meson propagator of the V-RMF and E-RMF models 
in stabilizing the transversal modes of the SNM for q ~ 0.5 fm and g ~ 5 fm can be seen in 
Fig.H 

If we extend this analysis to the asymmetric nuclear matter (ASM) case (Fig. [5]), it is 
found that different with the unstable regime of transversal mode, the instability of lon- 
gitudinal mode depends sensitively on the proton fraction in nuclear matter and after a 
certain critical proton fraction, which is less than the neutron fraction, this unstable regime 
disappears. 

Thus the longitudinal mode has a strong correlation with the isovector sector of the model. 
In the case of multi-component matter (p, n, e, and p) in (3 stability condition, a similar 
instability trend of the longitudinal and transversal modes with the ones obtained in the 
SNM case has been met. This means that the nucleons (protons and neutrons) play the main 
role for determination of the instability. However, a fine tuning in the isovector sector 
can shift the instability of longitudinal mode to higher critical density and momentum. 
Details of these results will be reported elsewhere. 

In conclusion, we have studied nuclear matter instability caused by particle-hole excita- 
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tions at go = in the RMF models at high density for longitudinal and transversal modes. 
It is found that in both modes the unstable regimes are very sensitive to the parameter 
set used. This opens the possibility to use the instability analysis as a tool to explore the 
applicability of parameter sets of the existing RMF models in high density applications. For 
example, in transversal modes, without the presence of vector and/or vector-scalar nonlin- 
ear terms with particular parameter values, the unstable regimes can not vanish. In general 
additional nonlinear terms in the form of vector and vector-scalar coupling terms of the V- 
RMF and E-RMF models improve the stability of both models at high densities. Finally we 
have observed that the longitudinal mode is sensitive to the isovector sector of the model. 
Thus, some important obscrvables for neutron stars, like proton fraction and asymmetry 
energy, can be related to this mode. 
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